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The average time T required for high-fidelity readout of quantum states can be significantly 
reduced via a real-time adaptive decision rule. An adaptive decision rule stops the readout as soon 
as a desired level of conhdence has been achieved, as opposed to setting a fixed readout time tf. 

The performance of the adaptive decision is characterized by the “adaptive-decision speedup,” tf fT. 

In this work, we reformulate this readout problem in terms of the first-passage time of a particle 
undergoing stochastic motion. This formalism allows us to theoretically establish the maximum 
achievable adaptive-decision speedups for several physical two-state readout implementations. We 
show that for two common readout schemes (the Gaussian latching readout and a readout relying on 
state-dependent decay), the speedup is bounded by 4 and 2, respectively, in the limit of high single¬ 
shot readout hdelity. We experimentally study the achievable speedup in a real-world scenario by 
applying the adaptive decision rule to a readout of the nitrogen-vacancy-center (NV-center) charge 
state. We find a speedup of « 2 with our experimental parameters. In addition, we propose a simple 
readout scheme for which the speedup can, in principle, be increased without bound as the fidelity 
is increased. Our results should lead to immediate improvements in nanoscale magnetometry based 
on spin-to-charge conversion of the NV-center spin, and provide a theoretical framework for further 
optimization of the bandwidth of quantum measurements. 
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I. INTRODUCTION 

Efficient discrimination of quantum states is impor¬ 
tant for, e.g., fast qubit readout rapid feedback 

and steering IMZ], preparation of nonclassical states of 
light [SJ IBlIlOj. and nanoscale magnetometry mi [n]. 
Given sufficient information about the statistics and dy¬ 
namics of a physical readout apparatus, it is possible to 
speed up a readout through a real-time adaptive decision 
rule (described below) [TO]. An adaptive decision 
rule can result in a significant reduction in the average 
time per measurement without any significant deteriora¬ 
tion of the readout fidelity. Adaptive decisions can be 
used to improve any readout scheme, in principle, given 
the ability to continuously monitor the state on a time 
scale that is short compared to the time required for high- 
fidelity readout [Tl|21|^ . 

Adaptive-decision problems have a long history in 
probability theory. Their mathematical root can be 
traced back as far as the famous “Gambler’s ruin” prob¬ 
lem introduced by Pascal in the I7th century [23]. A 
general mathematical theory of adaptive decision rules, 
known as sequential analysis, was developed during the 
World War II [53] and has since become a well-established 
part of statistical decision theory [25]. Implementing an 
adaptive decision rule requires the ability to update a 
stochastically varying measure of confidence in the state, 
typically a likelihood function, in real time. A mea¬ 
surement outcome is then chosen when the confidence 
measure first achieves a desired value. Adaptive-decision 
problems are thus closely related to first-passage time 
analysis [TO] Hz]; the stochastically varying confidence 
measure can be treated as the coordinate of a diffus¬ 
ing particle, which crosses a boundary when the desired 


confidence is reached. An adaptive decision allows for 
a reduction in the average measurement time, associated 
with a corresponding “adaptive-decision speedup.” Here, 
we formulate the adaptive decision for a two-state read¬ 
out in terms of a first-passage time problem. We use 
this formalism to theoretically establish the maximum 
achievable speedups for several physical readout models. 
Moreover, we demonstrate that significant speedups can 
be achieved in practice by experimentally characterizing 
the adaptive-decision speedup for the detection of a NV- 
center charge state, effectively doubling the bandwidth 
of such a measurement. 

The maximum achievable speedup is ultimately deter¬ 
mined by the specific dynamics of a given readout appa¬ 
ratus. Some readouts may allow for a very large speedup, 
while others are fundamentally limited. Here, we obtain 
upper bounds for the adaptive-decision speedup of two 
commonly encountered readout schemes. The first model 
we consider is the paradigmatic Gaussian latching read¬ 
out, in which the two states are distinguished by two 
noisy signals of constant but distinct intensity [28l [29] . 
This model can be a good approximation of, e.g., a read¬ 
out for the singlet and triplet states of two electron spins 
in a double quantum dot m, readouts of electron spin 
states based on repetitive measurement of nuclear spin 
ancillae [T71[TO], or a readout of superconducting qubits 
coupled to a microwave cavity [55] . The second model we 
consider is the readout based on state-dependent decay, 
in which the state is identified conditioned on a measured 
decay event. Readouts of, e.g., trapped ion qubits [3], 
NV-center spin qubits mm, and semiconductor spin 
qubits [T3j [16] can be approximated by this model. In 
the case of the Gaussian latching readout, we find that 
the speedup is bounded by a factor of 4 as the fidelity 
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is increased, while for the case of state-dependent decay, 
the speedup is bounded by a factor of 2. 

To experimentally study the achievable speedup, we 
implement an adaptive decision rule for a readout that 
discriminates between two charge states of a single NV 
center in diamond [121 EQ] . The NV-center charge read¬ 
out relies on the ability to distinguish strong from weak 
fluorescence upon illumination of the NV-center impu¬ 
rity. Remarkably, the NV-center charge readout can ap¬ 
proach either of the readout models described above in 
distinct limits. However, for realistic experimental pa¬ 
rameters, the dynamics of the NV-center charge readout 
will be intermediate between these two extremes. We 
therefore analyze and experimentally quantify the asso¬ 
ciated speedup. We update our level of confidence in 
the charge state using a quantum trajectory formalism 
that can easily be generalized to more complex read¬ 
out schemes. We find an adaptive-decision speedup « 2 
both for our experimental parameters and the experi¬ 
mental parameters of Ref. |12j , in which a similar charge- 
state measurement has recently been performed. Since a 
shorter average measurement time increases the number 
of measurements that can be performed on a NV-center 
spin in a given time, this result should directly improve 
the sensitivity of magnetometry based on spin-to-charge 
conversion of the NV-center spin m- A similar approach 
can, in principle, be applied to magnetometry using the 
standard spin readout of the NV center (SHHIS]- 


An adaptive decision can only speed up the aforemen¬ 
tioned readout schemes by a factor of order unity. How¬ 
ever, we show that the speedup can become parametri¬ 
cally large for other schemes. More precisely, we propose 
a simple readout, based on the discrimination of two dis¬ 
tinct decay channels, for which the speedup becomes un¬ 
bounded as the fidelity is increased. Such a readout could 
be realized in either atomic or quantum-dot systems. 


This article is organized as follows. In Sec. jlTj we de¬ 
scribe the general features of an adaptive decision rule 
that uses the likelihood ratio as a measure of confidence. 
In Sec. |HI[ we discuss the ma ximal speedup for the Gaus¬ 
sian latching readout (Sec. HI A) and for the readout 
based on state-dependent decay (Sec. HIB). We show 
that the speedup is bounded by a factor of 4 for the Gaus¬ 
sian latching readout and by a factor of 2 for the state- 
dependent decay readout. In Sec. TO we introduce the 
charge readout of the NV center (Sec. IV A). We describe 
our experimental setup as well as how the parameters 
of the charge dynamics were extracted from experimen¬ 
tal data (Sec. jlV B ). We then discuss how to apply the 
adaptive decision rule to the NV-center charge readout 
(Sec. [!v^ . We assess the performance of the adaptive 
decision rule both by performing Monte Carlo simula¬ 
tions based on experimental parameters and by applying 
the adaptive decision rule directly to experimental data. 
We demonstrate a factor « 2 speedup for the readout of 
the NV-center charge state. In Sec.jV) we discuss a read¬ 
out scheme where the speedup increases without bound 
as the readout fidelity is increased. We summarize in 


Sec. VI Technical details are given in Appendixes P|E| 


II. ADAPTIVE DECISION RULE 

The goal of readout is to discriminate between two 
particular states of a system of interest, |-|-) and |—). To 
achieve this, we typically let the system interact with a 
measurement device for a finite amount of time t, the 
readout time. During the time t, the measurement de¬ 
vice records a time-resolved trajectory in the form of, 
e.g., an electrical or photonic signal. As illustrated in 
Fig-Qa), a distinct trajectory is recorded when the state 
is either |-|-) (blue) or |—) (red), making it possible to 
discriminate between the two states. In the presence of 
noise, perfect discrimination is not possible. We then de¬ 
cide whether the state was most likely |-|-) or |—) based on 
the entire trajectory acquired during the readout time. 
For an optimal readout, this decision maximizes the read¬ 
out fidelity F, namely, the probability that the state is 
identified correctly. 

There is a tradeoff between readout fidelity and read¬ 
out time; a larger fidelity is typically achieved at the cost 
of a longer readout time. To achieve a given fidelity tar¬ 
get, the simplest approach is to choose a fixed readout 
time tf from measurement to measurement. The time 
tf is then selected to achieve the desired fidelity. An al¬ 
ternative strategy is to implement an adaptive decision 
rule. In this approach, each individual measurement is 
stopped as soon as a stopping condition is satisfied. The 
stopping condition is associated with reaching a mini¬ 
mum level of confidence in the state. A typical example 
of a stopping condition is illustrated in Fig. [^b). When 
the fidelity is high, a large subset of the possible trajec¬ 
tories will achieve a high level of confidence in a short 
amount of time, while others will require substantially 
longer times. Therefore, the average time T required for 
high-fidelity adaptive readout of the state may be signif¬ 
icantly shorter than the fixed time tf required to achieve 
the same fidelity with a nonadaptive readout. This is 
the main idea underlying the field of sequential detection 
theory, which aims to discriminate between competing 
hypotheses with the shortest possible average sampling 
time [Ml ESI 1311 [33]. The relative performance of the 
adaptive decision rule is characterized by the adaptive- 
decision speedup, tf/T. In this work, we analyze the 
maximum achievable speedup for three different read¬ 
out schemes listed in Table [H Each scheme will be de¬ 
scribed in detail in the following sections. Note that the 
adaptive decision rule discussed in this work does not 
require measurement feedback. However, it can be com¬ 
bined with measurement feedback to further reduce the 
average readout time m- 

We now formalize the ideas introduced above. For a 
given readout time t, the noisy state-dependent trajec¬ 
tory can be represented by a function defined in 

the interval t' G [0, t], as depicted in Fig. Ba). Given the 
observed trajectory, we wish to decide whether the initial 
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Readout scheme 


Maximal t//T 


Gaussian latching readout 
State-dependent decay 
Decay-channel discrimination 


4 

2 

oo 


TABLE I. Summary of the maximal speedups for the three 
readout schemes discussed in this work. 


state was most likely |-|-) or |—) in a way that maximizes 
the fidelity F [36j . Alternatively, we wish to minimize 
the error rate e = 1 — F. To do so, we calculate the like¬ 
lihood ratio At = P(i/)t|-l-)/F(V't|—), or equivalently the 
log-likelihood ratio At = In At. Here, P(^t|±) is the prob¬ 
ability of observing the trajectory ^t given that the state 
is initially |±). The log-likelihood ratio expresses a level 
of confidence in the state. When the prior probability 
for |±) is P(±), the decision rule that minimizes the er¬ 
ror rate is to choose the maximum a posteriori estimate 
(MAP) of the state; if At > Ath (At < Ath), we decide 
that the initial state was most likely |-|-) (|—))• Here, 
Ath = In [P(—)/P(+)] is the optimal decision threshold. 
When the prior probabilities are assumed to be equal, 
P{+) = P{—) and Ath = 0, the decision reduces to the 
maximum-likelihood estimate (MLE). In most of this ar¬ 
ticle, we assume equal prior probabilities so that the MLE 
is optimal. We briefly discuss unequal prior probabilities 
and the MAP in Sec. IIV C 31 

When using a nonadaptive decision rule, we fix t = ty 
for each measurement and base our decision on the value 
of Xtf (nonadaptive MLE or MAP). When using an adap¬ 
tive decision rule, in contrast, we stop each measurement 
as soon as At > A+, At < A_ or t > tM, as illustrated 
in Fig. Bb). Here, A± are stopping thresholds and tjvf 
is a maximum readout time. We then base our decision 
on the value of At at the stopping time (adaptive MLE 
or MAP). The average value of t from measurement to 
measurement gives the average readout time T. For each 
of the three readout schemes listed in Table [l] , we have 
proven that a symmetric choice of stopping thresholds, 
A_|_ = —A_ = A, minimizes T for a given error rate e when 
the states are equally likely. For the readout discussed in 
Sec. |IV[ we optimize the stopping thresholds numerically. 



Time, t' 



Readout time, t 


FIG. 1. (a) Readout trajectories for the Gaussian latching 
readout. The trajectory for state |-|-) (solid blue) has an av¬ 
erage of -1-1, while the trajectory for state |—) (dashed red) 
has an average of —1. Both trajectories have the same signal- 
to-noise ratio per unit time r. The identification of the state 
is made based on the entire trajectory acquired during the 
readout time t (dotted black), (b) Schematic representation 
of the adaptive decision rule for the Gaussian latching read¬ 
out. As soon as the log-likelihood ratio At (solid magenta) 
satisfies the stopping condition At > A, At < —A or t > tm, 
data acquisition is stopped. The sign of At is then used to 
choose the qubit state (dashed blue for |-|-) and dotted red 
for |—)). Here, we choose |-|-). 


A. Gaussian latching readout 


III. BOUNDED SPEEDUP 


In this section, we apply the ideas of Sec.|TT]to two com¬ 
monly encountered readout schemes, the Gaussian latch¬ 
ing readout and the readout based on state-dependent 
decay. Although these schemes are usually not perfectly 
realized in practice, they are relevant to understand the 
limitations of a wide variety of readouts. In particular, 
both models can be mapped to extreme limits of the NV- 


center charge readout discussed in Sec. IV 


The Gaussian latching readout is the simplest, most 
widespread and most tractable model of state discrim¬ 
ination [2H]- In this scheme, each state |±) gives 
rise to a state-dependent trajectory with constant av¬ 
erage ±1 and subject to Gaussian white noise, as il¬ 
lustrated in Fig. Ba). A latching readout is charac¬ 
terized by the absence of state relaxation at all times. 
Even though most readouts are limited by state relax¬ 
ation [saisTi I38j . many implementations will approxi¬ 
mately become latching readouts as the error rate is de¬ 
creased 

Formally, the probability distributions for ipt{t') con- 
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ditioned on the state are 

P(V't|±) = Aexp|-^y dt'[V't(t') T 1]^|- (1) 

Here, r is the power signal-to-noise ratio per unit time 
and H is a normalization constant independent of the 
state. From Eq. 0 . it is straightforward to show that 
the log-likelihood ratio is At = 2r J* dt' 'tpt{t'). Thus, in 
the case of the Gaussian latching readout, a maximum- 
likelihood decision can be made by integrating the tra¬ 
jectory up to time t. 

If the trajectory is acquired over the fixed interval 
[0,t/], the error rate takes the form [29]: 



where erfc is the complementary error function |39j . A 
derivation of Eq. ([^ is included as part of Appendix 
Alternatively, we can implement an adaptive decision 
rule. Because of the symmetry of the problem under 
interchange of |-|-) and |—), it is optimal to choose sym¬ 
metric stopping thresholds. Therefore, we stop the read¬ 
out as soon as A* > A, A* < —A or t > tM- As in the 
nonadaptive case, we choose the state according to the 
sign of A( at the end of the sequence. The adaptive de¬ 
cision rule is illustrated in Eig. 0b). The error rate e 
and average readout time T are then defined paramet¬ 
rically as a function of both A and Im ■ The calculation 
of e and T can be recast as a first-passage time prob¬ 
lem [HIITI, as detailed in Appendix]^ For clarity, here 
we only give the result for the case when the probability 
of reaching t > is negligible, rtM ^ min (l, A^) (see 
Appendix]^. The corresponding error rate and average 
readout time are given by 


e 


1 

1 -I- ’ 


T = — tanh 
2r 



(3) 


By varying the stopping threshold A, we can map the 
functional relation e(T). As e —> 0, Eq. ([^ implies that 
e « exp(—2rT). In the nonadaptive case, Eq. ([^, we 
instead have e « exp(—r t//2)/ y/2-Kr tf. Therefore, in 
the limit of small e, the adaptive decision rule asymptot¬ 
ically achieves the same error rate at an average time 4 
times shorter than the nonadaptive decision rule [25] [34] . 
Thus, the maximal speedup is given by 

y~4, (e^O). (4) 

Here and throughout this article, indicates an 

asymptotic equality. 


B. State-dependent decay 

As a second example, we consider the class of readouts 
relying on state-dependent decay. In these schemes, the 



Readout time, t 


FIG. 2. Adaptive decision rule for the state-dependent decay 
readout. The readout is stopped as soon as the likelihood 
ratio At (solid magenta) satishes At < —A or At > A. In this 
case, this pair of conditions is equivalent to stopping the read¬ 
out as soon as an event is detected before a maximum time 
tM = At (dot-dashed black). We choose the state according 
to the sign of At at the stopping time (dashed blue for |-|-) 
and dotted red for |—)). In this example, we choose |—). For 
the magenta trajectory shown above, this choice results in an 
error since an event occurs after time Im, as indicated by the 
sudden divergence in At. 


state |-|-) decays to |—) with probability per unit time 
T~^. The detection of a decay event, e.g., the detection 
of an emitted photon or of a tunneling electron, indicates 
that the initial state was |+). In contrast, the absence 
of a decay event indicates that the state was |—). This 
situation is approximately realized in, e.g., trapped ion 
qubits [1], NV-center spin qubits El [IB], and semicon¬ 
ductor spin qubits [THU]. To facilitate the discussion, 
we assume that the events are detected with perfect ef¬ 
ficiency. Moreover, we assume that a detected event is 
always the result of a transition; i.e., there are no “dark” 
counts. In this case, the log-likelihood ratio is At = —t/r 
when no event has been detected after a readout time 
t because the detection probability of a decay event de¬ 
creases exponentially with time. When an event is de¬ 
tected, At suddenly becomes inhnite because an event is 
guaranteed to have resulted from the initial state being 
|-|-). This scenario is illustrated in Fig. 

When using a nonadaptive decision rule, we wait for 
a time tf > t and base the decision on whether a decay 
event is detected in the interval [0,t/], as indicated by 
the sign of At. The error rate is given by the probability 
that the event occurs after time tf assuming an initial 
state |-|-), multiplied by the probability 1/2 that 

the state was initially |-|-): 



In contrast with the nonadaptive case, the adaptive 
decision rule stops the readout as soon as A* < —A or 
At > A, as shown in Fig.[^ Note that, for this particularly 
simple readout scheme, the use of a maximum readout 
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time tM is redundant since the above stopping condition 
is the same as stopping the readout as soon as an event 
is detected before a maximum time Im = A r. Variations 
of this adaptive decision rule have been implemented in 
Refs. [HIH]. The error rate is now given by the probability 
of an event occurring after time Im- To achieve the same 
error rate as the nonadaptive decision rule, Eq. , we 
thus set tM = tf. However, we find that the readout 
time is now given, on average, by (1 — when 

the state is |+). Thus, the average readout time is T = 
[(1 — + tf] /2 because, as noted in Ref. [4], we 

must always wait for the full duration tf when the state 
is |—). Thus, once again, we find a bounded speedup in 
the limit of vanishing error rate, t/ —> oo: 

(e^O). (6) 


IV. EXPERIMENT: NV-CENTER CHARGE 
READOUT 

The analysis of Sec. |III| provides limits on the achiev¬ 
able speedup for two idealized readout mechanisms. 
Many experimental realizations of a two-state readout 
can approach one or the other model in different limits, 
but under typical experimental conditions neither limit 
can be strictly realized. We must therefore use an algo¬ 
rithm for the adaptive decision rule suited to the under¬ 
lying readout mechanism. 

Here, we provide an explicit implementation of an 
adaptive decision rule for experimental two-state read¬ 
out in a regime that lies intermediately between the two 
idealized examples of Sec. |HI[ Specifically, we examine 
fluorescence-based charge-state detection of the NV cen¬ 
ter in diamond. As explained below, this system ap¬ 
proaches a Gaussian latching readout in the limit of high 
fluorescence rates and long charge-state relaxation times, 
while it can begin to resemble state-dependent decay for 
very low excitation powers. In the following, we describe 
the main features of the NV charge readout system and 
use experimental data to extract the underlying system 
parameters. We then both simulate and experimentally 
implement the adaptive decision rule with the help of a 
simple algorithm that uses our knowledge of the dynam¬ 
ics to update the likelihood ratio in real time. For our 
experimental parameters, we hnd a speedup oitf jT « 2. 
The resulting improvement in detection bandwidth can 
improve the sensitivity of an NV-center magnetometer 
using spin-to-charge conversion m- 


A. Charge dynamics of the NV center in diamond 

The NV center in diamond is typically observed in 
two charge states, namely, the negatively charged NV“ 
and the neutral NV®. In keeping with our notation, 
we label NV“ and NV® by |-|-) and |—), respectively. 


When the NV“ impurity is illuminated with yellow light 
(« 594 nm), transitions between the NV“ ground and 
excited states scatter photons which are detected at a 
rate 7 +. The frequency of the yellow light is below the 
minimum threshold required to resonantly excite optical 
transitions of NV®, and thus only a residual detection 
rate 7 _ ^ 7 + is observed in the neutral state. Note 
that the rates 7 ± include the effect of dark counts and 
imperfect detection efficiency. This difference in ffuo- 
rescence rates enables direct readout of the NV-center 
charge state muD]. In addition, the incident light can 
cause ionization (recombination) of the NV“ (NV®) state 
with probability per unit time r_|_ (r_). As illustrated in 
Fig-i a), this leads to an alternating process of ioniza¬ 
tion and recombination that limits readout fidelity [40| . 
Typically, the condition min( 7 _|_, 7 _) ^ max(r_|_,r_) is 
satisfied, so many transitions occur between ionization 
and recombination events. The trajectory ij^t can there¬ 
fore be modeled as a hidden two-state Markov process 
subject to state-dependent Poissonian noise of intensity 
7 ± [121 [20]. Because of detector bandwidth limitations, 
the maximum readout time tM is usually separated into 
N bins of duration St = tM/N, so we observe Srii pho¬ 
tons in bins i = 0,1, 2, ... , V — 1. Thus, the measured 
trajectory can be represented as an V-component vector 
t/’jv = (<5nAr_i, ... ,Sno). 

Note that the charge readout of the NV center ap¬ 
proaches the two cases discussed in Sec. HI in different 
limits. When min( 7 _|_, 7 _) ^ St~^, the noise statistics of 
the NV-center charge readout are approximately Gaus¬ 
sian, albeit with asymmetric signal-to-noise ratios per 
unit time: 


r± = 


( 7 + - 7)2 


47 ± 


(7) 


To perform a high-fidelity readout, the rate of in¬ 
formation gain must be faster than state relaxation, 
min (r+, r_) ^ max (r+ , r_ ). For an adaptive advan¬ 
tage, the rate of information gain must also be slower 
than the sampling rate, max(r_|_,r_) St~^. Combined 

with min( 7 +, 7 _) 3> the last inequality implies 

the necessary condition max(r_|_,r_) ^ min( 7 _|_, 7 _) or 
I 7 + — 7 _| ^ 2 min( 7 +, 7 _). We conclude that when 
state relaxation is negligible and when 7 _|_ and 7 _ are 
both large and comparable in magnitude, the NV-center 
charge readout is approximately described by the model 
of Sec. HI A with signal-to-noise ratio r « r+ « r_. In 


contrast, when 7 + ^ max( 7 _,r+,r_), the NV-center 
charge readout can be appro ximate ly modeled by the 


state-dependent decay of Sec. 


HIB 


with T = 


1 + 


In¬ 


deed, in this limit the presence or absence of a single pho¬ 
ton on a time scale of a few 7 )^^ is sufficient to choose the 
state with high confidence. Therefore, the detection of 
additional photons adds little information, and the state- 
dependent decay model is a good approximation. For our 
parameters, the experiments described below do not fall 
strictly into either of these limits. It is therefore inter¬ 
esting to characterize the speedup in this intermediate 
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(c) 


594 nm 


30 seconds 


532 nm 
Counting 


50 ms 1^ 

JirUlR - JlflflfUTrackingL 

-►i i-*- step 

100 jis 


the time scale for ionization and recombination, and each 
time bin (100/xs) is much shorter. A subset of such a 
time trace is illustrated in Fig. [^a) with data rebinned 
to 10 ms. The experiment is repeated to obtain statis¬ 
tics on fluorescence and ionization rates, interleaved with 
a tracking step (532-nm green laser) that ensures the 
sample does not drift relative to the focal point of the 
microscope. The experimental sequence is illustrated in 
Fig-d c). Such data were acquired for excitation powers 
ranging from 0.55/xW to 5^W (as measured going into 
the objective). For comparison, the saturation intensity 
at 594 nm is estimated to be 2.5 mW. We verified that 7 + 
and 7 _ scale as the laser intensity and that F-i- and r_ 
scale as the laser intensity squared, as shown in Ref. [20] . 
Data sets were measured for two (111) oriented NV cen¬ 
ters, which yielded similar results. Very weak excitation 
is used to ensure that 7 + and 7 _ are much larger than 
F+ and F_. 

For definiteness, we restrict the following analysis to 
a single power setting. Specifically, the data sets for the 
lowest laser power setting of 0.55 fiW were used to extract 
the parameters 7 ± and r±, as detailed in Appendix |E 1[ 
We obtain the rates 7 + « 720 Hz, 7 _ « 50 Hz, r+ « 
3.6 Hz and r_ « 0.98 Hz, which are required to construct 
an algorithm for an adaptive decision rule. 


FIG. 3. Schematic representation of (a) the NV-center charge 
dynamics described in Sec. IV A[ (b) the experimental appa¬ 
ratus and (c) the experimental sequence. 


regime. 


B. Experimental setup and data 

We characterize the performance of the adaptive de¬ 
cision rule using experimental data. The data were 
acquired using a home-built confocal microscope with 
594-nm excitation. The 594-nm wavelength lies in be¬ 
tween the zero-phonon line of the NV“ (637 nm) and 
NV° (575nm), so it only efficiently excites the negatively 
charged state. The excitation light is focused through a 
high numerical aperture objective (NA 1.35) onto a (111) 
cut chemical-vapor-deposition-grown diamond in which 
single defects can be resolved. Photons emitted in the 
wavelength range 645 — 800 nm are collected and detected 
with a single-photon counter. This detection range over¬ 
laps strongly with the NV“ fluorescence spectrum and 
only weakly with the NV° fluorescence, while rejecting 
Raman scattering from the diamond. The experimental 
setup is illustrated schematically in Fig. |^b). 

To acquire fluorescence trajectories, 594-nm excitation 
is applied continuously and photon counts are recorded 
for 30 s in bins of 100 /xs, with a small separation of 8.3 ns 
necessary for our field-programmable gate array (FPGA) 
card to process the counts. These time scales are chosen 
such that the total duration (30 s) is much greater than 


C. Readout error analysis 

In Refs, [miini, single-shot readout of the charge state 
was performed by counting the total number of photons 
detected in the interval [0,t/] (see Appendix [Pj). In con¬ 
trast, implementing the adaptive decision rule requires 
an efficient algorithm to update the likelihood ratio after 
each time bin [IIMIIIZIIIIHIII- A hidden-Markov-model 
algorithm suitable to the NV-center charge readout has 
been developed in, e.g.. Refs. |151|33j. Here, we use an 
equivalent quantum trajectory approach, which is more 
easily generalized to quantum systems with multiple lev¬ 
els and coherent internal dynamics. Moreover, we find 
an update rule that is valid for an arbitrary value of the 
bin size 6t. In particular, more than one ionization or 
recombination event may occur during one time bin, an 
occurrence that becomes more likely at high illumination 
power |5D| . The update rule of Ref. jUj is recovered when 

r±st < 1. 

As shown in Appendix the likelihood ratio Ajv = 
P('0Ar|+)/P(V'Ar|—) after N bins is given by 


Tr 

1 

. M((5no)4+^ 

Tr 

M((5nAr_i) . 

. M(fco)4”^ 


In Eq. ([^, = (1,0)^ and = (0,1)^ are ini¬ 

tial state vectors representing |-|-) and |—), respectively. 
Here, the trace Tr of a vector is defined as the sum of its 
elements. The matrix M(5 ?t,) is a 2 x 2 update matrix for 
the detection of Sn photons in a bin. The exact form of 
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FIG. 4. (a) Experimental photon-count trajectory illustrating 
the state-dependent fluorescence rate and two-level switching 
of the NV-center charge state for an illumination power of 
0.55 ^W. Here, the photon detection rate fluctuates between 
7 + ~ 720 Hz for NV“ and 7 _ « 50 Hz for NV°, as discussed 
in Sec. |IV B| The corresponding average counts per 10-ms 
bin are indicated by the horizontal dotted black lines. The 
NV-center ionization and recombination rates are found to 
be r+ « 3.6 Hz and T- ~ 0.98 Hz, respectively, (b) Log- 
likelihood ratio At (solid magenta) and counts per 100-/is bin 
(dashed red) for the first 60 ms of the data shown in (a). The 
log-likelihood ratio was calculated using Eq. ^ and the ex¬ 
perimental parameters given in Sec. |IVBl The discontinuities 
in At indicate the detection of a photon. 


the update matrices is given in Appendix]^ They can be 
stored for all Sn < Sumux, where is the maximum 

number of photons with a non-negligible probability to 
occur. The likelihood ratio, Eq. can then be updated 
in real time through simple matrix multiplication. An ex¬ 
ample of the application of Eq. ([^ to experimental data 
is given in Fig. Sb). Note that IVE is nothing but the 
measurement superoperator for a classical Markov pro¬ 
cess subject to Poissonian noise. The approach can thus 
be generalized to any Markovian quantum trajectory by 
substituting M with the appropriate measurement super¬ 
operator for direct detection and the vectors by any 
pair of initial states to be discriminated. In particular, 
the standard spin detection of the NV center modeled in, 
e.g., Refs. pSl - HT] can be processed in this way. 


In the following, we first use Monte Carlo simulations 
to assess the ideal theoretical performance of the adaptive 
decision rule. We then apply the adaptive decision rule 
to experimental data to demonstrate a speedup under 
real experimental conditions. In all cases, we choose the 
experimental time bin St = 0.1ms and a maximum ac¬ 
quisition time tjvf = = 25 ms. We use this value of St 

and the rates given in Sec. |IVB] to calculate and store the 
update matrices M((5n) for Sn = 0,1,... ,Sn^ax- Here, 
Sn^ax = 5 is chosen to match the maximum number of 
photons observed in any one time bin of the experimen¬ 
tal data sets. The residual probability of observing more 
than Sn^ax photons in one time bin is « 10“^°. 


1. Monte Carlo simulations 

For both initial states |±), we use the procedure de¬ 
scribed in Appendix [^to generate 1 x 10® random trajec¬ 
tories i/'Af in the interval [0, NSt]. We then apply both the 
nonadaptive and the adaptive MLE to each trajectory 
to obtain a Monte Carlo estimate of the average read¬ 
out time T and error rate e, as detailed in Appendix [C| 
The resulting error-rate curves are shown in Fig. For 
comparison, we also show the error rate of the photon¬ 
counting method modeled in Ref. [12j and summarized 
in Appendix]^ In this approach, the total number of de¬ 
tected photons is compared to a threshold to choose the 
state. As expected, the nonadaptive and adaptive MLE 
achieve a somewhat lower minimum error rate than pho¬ 
ton counting. More importantly, however, we see that 
the adaptive MLE achieves the minimum error rate of 
the photon-counting method in approximately half the 
time necessary for the nonadaptive MLE to achieve the 
same goal. Thus, the bandwidth of the already efficient 
NV-center charge readout is doubled through signal pro¬ 
cessing alone. We have verified that a similar advantage 
also exists for the experimental parameters of Ref. m- 
Since this speedup enables a larger number of NV-center 
charge measurements to be performed in a given amount 
of time, it should lead to a substantial improvement in 
the sensitivity of the spin-to-charge conversion magne¬ 
tometer in cases where the measurement duty cycle is 
limited by readout time. 

We note that the parameters given in Sec. |IVBl satisfy 
the condition for which the NV-center charge readout re¬ 
duces to the state-dependent decay readout of Sec. |IIIB 
T~^ = 7 - 1 - > max( 7 _,r+,r_) = 7 -. It is therefore plau¬ 
sible that the speedup « 2 obtained in Fig. has at 
least partially the same origin as that given in Eq. §. 
Indeed, the speedup « 2 observed in Ref. [3] was essen¬ 
tially explained by the state-dependent decay model for 
a readout whose dynamics are formally the same as the 
NV-center charge readout dynamics. However, the model 
of Sec. |HI 5] fails when the maximum readout time tM be¬ 
comes comparable to ~ 20 ms, resulting in a much 
higher saturation error rate in Fig. than that predicted 
by Eq. ([^ with r = « 1.4ms and t/ = = 25 ms. 



































FIG. 5. Error rate as a function of average readout time 
for the experimental parameters of Sec. |IV B[ The lines give 
the Monte Carlo simulated error rate, and the symbols give 
the experimental values corrected for preparation errors using 
Eq. (E71. We display the error rate for the photon-counting 
model discussed in Appendix P] (solid blue and blue circles), 
for the nonadaptive MLE (dashed orange and orange squares), 
and for the adaptive MLE (dotted green and green triangles). 
The horizontal dotted black lines indicate the minimum error 
rate of the photon-counting (1.9%) and MLE methods (1.5%). 
As indicated by the vertical dotted black lines, the speedup 
of the adaptive MLE compared to the nonadaptive MLE is 
tfjT « 1.9 when the target error rate is the minimum error 
rate of the photon-counting method. 


In our experiment, it is therefore necessary to use the 
update rule of Eq. (|^ to obtain an accurate estimate of 
the error rate. In a context where it may be favorable or 
necessary to move towards the regime of the high-fidelity 
Gaussian latching readout of Sec. |III A[ speedups larger 
than 2 should be possible. Entering this regime may not 
be advantageous for the NV-center charge readout since 
the increase in excitation intensity required to reach the 
Gaussian regime results in a detrimental increase in the 
ionization and recombination rates. However, such an 
advantage may be possible in other systems with similar 
dynamics. 


parameters and Eq. (|^ are then used on the testing set to 
verify the model. The details are given in Appendix |E 2| 
We first perform a preliminary verification by com¬ 
paring the experimental distribution of log-likelihood ra¬ 
tios to the distributi on p redicted by Monte Carlo simula¬ 
tions [see Appendix E 2 . We find a close agreement be¬ 
tween experiment and theory, showing that the model of 
charge dynamics discussed in Sec. IV A provides a good 
description of the statistics of the experimental trajec¬ 
tories. We can therefore use our model to perform ap¬ 
proximate preparation of the charge state in postselec¬ 
tion. We use this preparation to verify that the adaptive- 
decision speedup exists when the adaptive decision rule 
is applied to experimental data. From the testing set, 
we prepare 8307 (30693) trajectories of 25 ms with ini¬ 
tial states j-b) ([—)). We then read out the state for 
each trajectory using the photon-counting, nonadaptive 
MLE, and adaptive MLE methods (i.e. assuming equal 
prior probabilities). Comparing the result of the readout 
to the preparation allows us to estimate the experimental 
error rate e = (e+ -I- e-)/2, where e± are the error rates 
conditioned on the preparation in state |±). To better 
compare the experimental results with theory, we then 
fit the experimental error-rate curves to the theoretical 
prediction, acc oun ting for an additional preparation er¬ 
ror r] [see Eq. (E7|]. This procedure gives r] = 2.22% as 
a single fit parameter. We emphasize that this adjust¬ 
ment is a uniform transformation of the error rate for all 
times and hence does not affect the measured speedup. 
The experimental error rate e obtained after compensat¬ 
ing for preparation error is shown along with the error 
rate determined by Monte Carlo simulation in Fig. [the 
uncompensated error rate e is shown in Appendix |E 2| 
for comparison]. The experimentally measured error rate 
thus fits the theoretical prediction very well. Therefore, 
we conclude that the adaptive-decision speedup « 2 dis¬ 
cussed in Sec. |IV CT] persists when the adaptive decision 
rule is applied to real experimental data in spite of pos¬ 
sible systematic errors and imperfect modeling. 


3. Unequal prior probabilities 


2. Experimental verification 


In the simulations of Sec. IV C1 the statistics of 


the simulated fluorescence trajectories are perfectly de¬ 
scribed by the two-state model assumed for readout. In 
real experimental conditions, however, the charge dy¬ 
namics may deviate from this model. In this section, 
we show that the adaptive-decision speedup is robust to 
imperfections in our modeling by applying the adaptive 
decision rule to experimental data. To independently ver¬ 
ify our model of charge dynamics, we split our data sets 
into a calibration set and a testing set. The calibration 
set is used solely to extract the experimental values of 


the rates 7 ± and r± given in Sec. IV B The extracted 


A balanced probability distribution for the initial state, 
P{+) = P{—), is often desirable because it enables the 
extraction of one bit of information per measurement. 
However, in applications such as the magnetometry pro¬ 
tocol of Ref. [12] , the probability distribution of the ini¬ 
tial state is, in reality, unbalanced, P{+) yf P{—)- In gen¬ 
eral, the dependence of the adaptive-decision speedup on 
the prior probabilities T’(±) is nontrivial and may depend 
on the particularities of the readout dynamics. A general 
analysis lies beyond the scope of this work. Here, we nev¬ 
ertheless discuss how to account for unequal prior prob¬ 
abilities. Moreover, we show that the adaptive-decision 
speedup of the experimentally relevant NV-center charge 
readout discussed in Sec. |IV CTj persists even for unequal 
prior probabilities. We distinguish the case where P(±) 
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are unknown from the case where P{i) are known. 

When the prior probabilities are unknown, the best 
strategy is to calibrate the readout thresholds Ath and 
A± by assuming that the prior probabilities are equal. 
This leads to the MLE readout discussed in Secs. IIV C 11 
and |IV C 2[ Under this assumption, we obtain the er¬ 
ror rates and average times conditioned on the 
state |±). When the MLE readout is used on an unbal¬ 
anced sample of charge states, the error rate and average 
readout times are then given by e = P(-|-)e*^“''^ -|-P(— 
and T = P(-|-)r^+^ + respectively. The sim¬ 

ulated error rate as a function of average readout time 
in this scenario is shown in Fig. [^a) for the experimen¬ 
tal parameters of Sec. [rV^ and P{+) = 0.25. We also 
show the corresponding error-rate curve for the photon¬ 
counting method (see Appendix [d|. Using the minimum 
error rate of the photon-counting method as a reference, 
we find that a substantial speedup tf/T « 1.6 persists 
even when the prior probabilities are unknown. 

When the prior probabilities are known, they can 
be used to improve the MLE readout by adjusting the 
thresholds to directly minimize e = P(-|-)e^'''^ -|-P(— 
for constant T = P{+)T^P^ + P{—)T^~\ This leads to 
the MAP readout mentioned in Sec. [u] (see Appendix [C| 
for details). We note that in a given application, the 
prior probabilities can easily be determined by using the 
MLE to estimate the relative proportion of |-|-) and |—). 
The simulated error rate as a function of average read¬ 
out time when using the MAP is shown in Fig. [^b) for 
the same experimental parameters and P(-l-) = 0.25. We 
find an adaptive-decision speedup tf/T « 1.8 using the 
minimum error rate of the photon-counting method as a 
reference, similar to the value obtained in Sec. IV C 1 for 
P(-l-) = 0.5. Therefore, we conclude that the adaptive- 
decision speedup persists for significant deviations from 
a balanced initial-state distribution. 


V. PARAMETRIC IMPROVEMENT IN 
SPEEDUP 


From the examples of Secs. |III| and |IV[ we might be 
tempted to conclude that the speedup is fundamentally 
bounded by a constant of order unity. Here, we show 
that there exist readout schemes where the speedup can 
become arbitrarily large in the limit of low error rate. 

To show this, we consider a simple variation of the 
state-dependent decay readout analyzed in Sec. |IIIB] We 
now suppose that the states |-|-) and |—) both decay to a 
third state |0) through different decay channels R and L, 
respectively. The two channels may be associated with, 
e.g., different polarizations of emitted photons or two dif¬ 
ferent leads into which the electron of a double-quantum- 
dot charge qubit can tunnel (see Fig. [^. If the two de¬ 
cay channels can be discriminated, e.g., with the help of 
a polarization analyzer, the state |±) can be read out. 
For simplicity we assume that both states decay with the 
same probability per unit time 




Average readout time, T (ms) 

FIG. 6. Monte Carlo simulated error rate for an unbalanced 
initial charge distribution with P{+) — 0.25 when (a) the 
value of P{+) is unknown (MLE) and (b) the value of P{+) 
is known (MAP). In both cases, we show the error rate for 
the photon-counting method of Ref. [12] (solid blue), the non- 
adaptive readout (dashed orange) and the adaptive readout 
(dotted green). The horizontal dotted black lines indicate the 
minimum error rates of the photon-counting method (1.6% 
in both cases) and of the optimal methods (1.4% for the 
MLE and 1.3% for the MAP). The dotted vertical lines show 
that significant adaptive-decision speedups of tf/T ~ 1.6 and 
tf/T ~ 1.8 are obtained when using the MLE and MAP, re¬ 
spectively. The discontinuous features in (a) occur because 
the MLE uses a nonoptimal discrimination threshold when 
P(±) / 0.5. 


When using a nonadaptive decision rule, we wait a 
time tf and base the decision on whether a i? or L 
decay event has been detected. If no event is detected, 
we choose the state at random with equal probability. 
Assuming, for simplicity, that the channels R and L can 
be perfectly discriminated, the error rate is given by the 
probability that no decay event has occurred up to time 
tf, multiplied by the probability 1/2 that the 

random decision fails. Thus, the error rate is still given 
by Eq. (§. 

The adaptive decision rule, in contrast, stops the read¬ 
out as soon as either a decay event is detected or a max¬ 
imum time Im is reached. Following the same argument 
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(a) 


|-) = |e,L> l+> = |e.i?> 



|0) = Iff) 



states, 1+) = |e,i?) and |—) = |e,L), both decay to the 
ground state |0) = \g) by emitting photons with right 
and left circular polarizations, respectively. This situ¬ 
ation is depicted in Fig. A polarization analyzer 

could then identify both states in a time « r, leading 
to the parametric improvement of Eq. ([^. In another 
example, suppose that we want to discriminate between 
an electron being in the rightmost (|-|-) = \R)) and left¬ 
most (|—) = |L)) islands of a double quantum dot. We 
imagine that, as shown in Fig.j^b), the potential barriers 
between the dots and leads are lowered at the beginning 
of the readout phase. An electron in the right (left) dot 
will then tunnel out into the right (left) lead and be de¬ 
tected, with the left and right detectors playing the role 
of the polarization analyzers in the previous example. 
In perfect analogy with Fig. [^a), stopping readout as 
soon as an electron is detected in either lead instead of 
waiting for a fixed time tf then leads to a parametric 
improvement in speedup, Eq. Note that while the 
assumption of perfect detection efficiency is difficult to 
realize for photon detection, it is usually not a limitation 
for the detection of electron charges in quantum dots. 


VI. CONCLUSION 


FIG. 7. (a) Readout scheme based on state-dependent polar¬ 
ization of an emitted photon. When the state is |-|-) = \e,Ft) 
(|—) = |e, L)) in the excited subspace, a right (left) circularly 
polarized photon is emitted. Both states decay to the ground 
state |0) = |(/) on a time scale r. (b) Similar scheme for read¬ 
out of the charge state of a double quantum dot. When the 
electron is in the right (left) island |i?) {\L}), it is detected 
after tunneling into the right (left) lead. Both states tunnel 
into the leads on a time scale r, leaving the double quantum 
dot in the empty state |0). 


as for the nonadaptive decision rule, the error rate is 
e = To achieve the same error rate as the non¬ 

adaptive decision rule, we must therefore choose = tf- 
However, we find that the average readout time is now 
T = (1 — As the error rate decreases, tf —)■ oo, 

the average readout time tends to a constant t so that 

(e^O). (9) 

Therefore, the speedup increases without bound as the 
error rate e is decreased. This parametric improvement 
in speedup is achieved by transferring all the information 
about the state to the channel degree of freedom, R or L. 
This means that the readout no longer relies on discrim¬ 
inating decay from the absence of decay. We can thus 
arrange the readout so that both |—) and |-|-) decay on a 
time scale r. Therefore, it is no longer necessary to wait 
for the entire duration tf when the state is |—). 

Such a scheme could be applied in a variety of sce¬ 
narios. For example, suppose that two atomic excited 


In summary, we have established the maximum achiev¬ 
able “adaptive-decision speedups” for several physical 
readout models in the high-fidelity limit. To achieve this, 
we formulated the adaptive-decision problem in terms of 
a first-passage time problem. We obtained bounds for 
the speedup of two commonly encountered readout mod¬ 
els. Specifically, we have shown that the adaptive de¬ 
cision rule can speed up the Gaussian latching readout 
(Sec. [Ill A ) by up to a factor of 4 in the limit of high 
readout fidelity, while it can s peed u p readouts relying 
on state-dependent decay (Sec. IIIBl by up to a factor 
of 2. To study the achievable speedup in a real-world 
scenario, we applied the adaptive decision rule to a read¬ 
out of the charge state of a NV center in diamond using 
a quantum trajectory formalism that incorporates the 
NV-center charge dynamics and accounts for experimen¬ 
tal imperfections such as dark counts and imperfect col¬ 
lection efficiency. Although the NV-center charge read¬ 
out reduces to the two aforementioned models in distinct 
limits, its dynamics are in an intermediate regime for 
typical experimental parameters. We have shown that 
a significant speedup can be achieved under these con¬ 
ditions. Specifically, we found a speedup « 2 both in 
our experiment and using the experimental parameters 
of Ref. [H]. Finally, we have proposed a readout scheme 
that leads to an unbounded speedup as the fidelity is in¬ 
creased, in stark contrast to the common readout models 
discussed previously. This scheme relies on discriminat¬ 
ing between two distinct decay channels instead of dis¬ 
criminating between decay and absence of decay. We 
have further provided avenues to realize such a scheme 
in atomic or quantum-dot systems. 
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Our results are already applicable to a wide variety 
of systems, and they provide a direction for the opti¬ 
mization of measurement bandwidth in experimentally 
relevant readouts of quantum states. In particular, we 
have shown that magnetometry based on the NV-center 
spin-to-charge-conversion readout can be improved with 
currently achievable experimental parameters. Moreover, 
our results show that apparent limitations to the adap¬ 
tive speedup can be overcome by careful redesign of the 
readout dynamics. In this work, our goal was to ana¬ 
lyze the fundamental limitations of the adaptive-decision 
speedup for extreme limits of several physical readout 
schemes. We expect direct extensions of the first-passage 
time formalism developed here to allow for the derivation 
of analytical speedup bounds in the presence of, e.g., un¬ 
wanted state relaxation, dark counts, or imperfect pho¬ 
ton collection efficiency. On a more fundamental level, 
our analysis provides the necessary framework to study 
how the addition of measurement feedback and coher¬ 
ent readout dynamics can modify the maximal speedups 
discussed here. Moreover, this formalism may allow for 
the direct characterization of the achievable speedups for 
parameter estimation, or multiple-state discrimination in 
general, subject to experimentally relevant noise. 
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Appendix A: Analytical treatment of the adaptive 
Gaussian latching readout 


In this section, we sketch the derivation of the error 
rate and average readout time for both the nonadaptive 
and the adaptive decision rule applied to the Gaussian 
latching readout. According to Eq. 0. the expectation 
and covariance of the trajectory iptit') conditioned on the 
state |±) are 


E[V^,(0|±] = ±I, 

Cov = r~^S{t' - t"), 

for t' < t and t” < t. Because the noise is Gaussian, 
all higher cumulants vanish. Gorrespondingly, the ex¬ 
pectation and variance of the log-likelihood ratio At = 
2r J* dt' are E(At|±) = ±2rt and Var(At|±) = 
Art. Therefore, At is a simple drift-diffusion process dis¬ 
tributed according to 


P(At|±) = G(±)(At,t) 

1 


y/Sirr t 


exp 


(At T 

8rt 


(A2) 


Here, is the Green’s function that solves the 

associated drift-diffusion equation for At. 


For a fixed readout time tf, the state is chosen accord¬ 
ing to the sign of At^. The error rates conditioned on the 
initial state being |±), are thus 


e(+) = 


e 


(-) = 


f dXtj, P(Xtj,\+), 

J —oo 


dXtf P{Xtf |—). 


(A3) 


By symmetry of the problem, 

Eq. (A2), we can then calculate the error rate 

[e(+)V7-)] /2 = e(±): 


Using 
e = 


e = 



(A4) 


To assess the performance of the adaptive decision rule. 


we must calculate both the error rate e = [' 


= U(+) 


-)] /2 


and the average readout time T = +T^ /2, 

where is the average readout time conditioned on 
the state being |±). Using the symmetry of the prob¬ 
lem again, we have e = and T = . We may 

therefore assume that the state is |-|-) without loss of 
generality. We now reformulate the calculation of e and 
T as a first-passage time problem ( an alternative 

formalism is given in Ref. [S^). More precisely, let e{X, t) 
and T{X,t) be the error rate and average readout time 
conditioned on knowing that the log-likelihood ratio is A 
at time t. We wish to obtain e = e(0, 0) and T = T(0,0). 
Beca use drift diffusion is a Markov process, we can use 
Eq. (A21 to condition the values of e and T on their 
possible values at time t + St: 


/ OO 

dx St)e{X + x,t + St), 

-OO 

/ OO 

dxG'^+\x,St)r{X + x,t + St). 

-OO 


(AS) 


We then expand e and T around a: = 0 and St = 0 on 
both sides and keep all terms of order St. In the limit 
St —>■ 0, Eqs. (AS) then reduce to Kolmogorov backward 
partial differential equations EillTl of the drift-diffusion 
type [IS]: 


1 de{X,t) de{X,t) d'^e{X,t) 

2r dt dX dX^ ’ 

I dr{x,t) dT{x,t) d^rix,t) 

2r dt dX dX^ 2r 

The adaptive decision rule is implemented by setting ap¬ 
propriate boundary conditions. Since we have assumed 
that the state is |-|-), an error occurs only when we stop 
with At < 0. We must therefore have £(±A,f) = ^(tA) 
and e{X,tM) = d{—X). Similarly, the rernaining read¬ 
out time after stopping must vanish, T(±A,t) = 0 and 
T{X,tM) = 0. Here, 6{x) is the Heaviside step function. 
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To solve Eqs. (A6) analytically, we first take Im —t oo. 
The corresponding solutions £oo(A) and 7^ (A) are inde¬ 
pendent of t and satisfy the ordinary differential equa¬ 
tions: 


5eoo(A) d'^e^{X) _ 
dX dX^ 

dTooiX) , d^TooiX) , 1 


(A7) 


dX 


dX^ 




subject to the boundary co nditi ons eoo(±A) = 0 (tA) and 
7^(±A) = 0. Solving Eqs. (A7) gives Eq. ([^: 


Coo — £oo(0) — 


1 -I- e-^ 


A 


Too = 7^(0) = — tanh ( - ) • 


(AS) 


7+ 7+ 



7- 7- 


FIG. 8. Unraveling of the NV-center charge readout dynam¬ 
ics. The state |-|-) (|—)) represents NV“ (NV'^). Photons 
are detected at state-dependent rates 7 ±. Each detection in¬ 
creases the number of observed photons n. The rates r± 
describe ionization or recombination processes that are not 
associated with the detection of a photon. 


To obtain the solution for finite tM, we write e{X,t) = 
£ao{X) + 'r]{X,t) andT(A,t) = Taaj X) -\- ({X, t). Substitut¬ 
ing these expressions in Eqs. ( |A6[ ), we find that rj and C 
satisfy the homogeneous equations 

_1 dviX,t) ^ d-q{X,t) d‘^r]{X,t) 

2r dt dX dX^ ’ 

1 dgx^t) _ dgx^t) d\{x,t) 

2r dt dX dX^ 


Appendix B: Derivation of the update matrices 

In this section, we derive the form of the update ma¬ 
trices M(n) using a number-resolved quantum trajec¬ 
tory formalism [49l [50]. The NV-center charge state 
at any given time t can be described by a state vector 
p = (p+, p_)^, where p± is the probability of finding the 
charge state |±). The state vector obeys the equation of 
motion of a Markovian two-level fluctuator: 


subject to the boundary conditions ri{±X,t) = 0, 

= S{-X) - eoo(A) and C(±A,t) = 0, C(A,tM) = 
—Too{X), respectively. For each of the Eqs. (A9), we de¬ 
compose the boundary condition at time tM in terms of 
the basis of right eigenfunctions of the drift-diffusion op¬ 
erator d\ + df. The eigenfunctions are chosen to vanish 
for A = ±A. We then propagate each component back¬ 
wards in time to find r?(0, 0) and C(OjO)- This gives the 
exact analytical expressions: 


e = eoo+ ^ 

m—0 

oo 

rT = rT^+^B^e-“-^‘" 

m—Q 


(AlO) 


Here, we define 


471^(771 -b 1/2)2 _|_ y2 ’ 

IbTTp cosh (|)(-l)"^(w -b 1/2) 
[47r2(m -b 1/2)2 _|_ ^2j 2 
1 ^ 27r2(m -b 1/2)^ 


p{t)^Cp{t). (Bl) 

The state at time t = 0 is p(0) = Pq. Here, C is the 
Lindblad superoperator of the two-level fluctuator (in the 
basis {|+)(+|,|-)(-|}): 



To analyze the photon emission statistics, we resolve (or 
unravel) the state vector in the photon number: 

OO 

p{t) = '^e{n,t), (B3) 

n—0 

where £(n, t) is the unnormalized state vector after mea¬ 
surement of n photons in the interval [0,t]. More pre¬ 
cisely, the probability of detecting n photons after time 
t is 


P{n,t) = Ti £{n,t), (B4) 

where the trace Tr is defined as the sum of all elements 
in the vector. Our goal is to find the measurement su¬ 
peroperator M.{n,t) so that P{n,t) can be expressed as 

P{n,t) = Tr [A4(n, t)pg] . (B5) 


Equation (AlO) reduces to Eq. (A8) when rtM >• 
min(l, A^). 


We must first obtain the equations of motion for £{n, t). 
The dynamical evolution of £(n, t) is illustrated in Fig.|^ 
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At each time step, there are two possible types of tran¬ 
sition. Either an ionization or recombination occurs (n 
is unchanged) or a photon is detected (n is increased by 
1). Thus, £{n,t) obeys the following set of coupled rate 
equations: 

£{n,t) = C£{n,t) — K [t{n,t) — t{n — l,t)]. (B6) 

The initial condition is ^(0,0) = £q = pp. Here, K en¬ 
codes the state-dependent photon detection: 




7+ 0 \ 

0 7-j' 


(B7) 


Note that summing Eq. (B6) over all n recovers Eq. (Bl) 
after applying Eq. (B3), as required by conservation of 
probability. The effect of dark counts and imperfect de¬ 
tection efficiency is simply to modify the rates 7 ±. 

Equation (B61 can be solved by introducing the char¬ 
acteristic function (or Fourier transform) |50j 


€(x,^) = ^£(n,^)e“^ 


= — / dxi{x,t)e 


(B8) 


where y is a counting field. Substituting Eq. (B8) into 
Eq. (B6), we have 


= '^(x)€(x,t). 


(B9) 


where 


n{x) = C - K + e^^K. (BIO) 

Solving Eq. ( |B9[ ) and inverting the characteristic function 
then gives 


£{n,t) = A4{n,t)£o, (BH) 


where the measurement superoperator Ad(n,t) is 


1 

M(n,t) = —I 

1 

M(n,t) = - — 
n\ dz^ 


(B12) 


In the last line, we expressed the result in terms of the 
z-transform variable, z = [49]. We see that the 

measurement superoperators are generated by the ma¬ 
trix function Q{z,t) = 

We can now calculate the probability of the trajectory 
ipN = (i5nAr-i, ... , Sno) given the initial state as 


P{i’N\£o) = Tr [M((5njv_i) ... M{Sno) £o ], (B13) 

where M(( 5n) = Ad.{Sn,St) is the desired update matrix. 
Equation ( |B13 ) can be used to compare the likelihood 
function for two arbitrary initial states in particular 
the two charge states = (1,0)”^ and £q~^ = (0,1)^. 


If we choose a single time bin, 6t = tf, M(n) generates 
the photon number distributions given in Ref. |12j (see 
Appendix |D|) . Moreover, expanding M(fc) for small St 
yields the continuous-time measurement superoperator 
for direct detection m- Note also that Eq. ( |B13| ) is an 
exact solution of the continuous-time filtering equations 
discussed in Ref. |12|. 


Appendix C: Details of Monte Carlo simulations 


In this section, we describe the Monte Carlo simula¬ 
tions of the error rates in more detail. For a given ini¬ 
tial state Pq, we generate a trajectory iteratively using 
Eq. (|B13| as follows: 


1. Knowing the state p^. in bin k, calculate the corre¬ 
sponding photon probability distribution Pk(ji) = 
Tr [M(5n)pj,] for Sn = 0,1,2,... ,Sn 

max ■ 

2. Sample a value Suk at random from the distribu¬ 
tions Pk{Sn). 

3. Find the postmeasurement state Pk+i = 
M{Snk)Pk/Pk{Snk). 

4. Go back to step 1 with initial state p^+i- 


We use this procedure to simulate 1 x 10® trajectories 
using driniax = 5, St = 0.1ms, Im = 25 ms, and the rate s 
extracted from experimental data (see Appendix E 11. 
Each trajectory is then processed using both the non- 
adaptive and adaptive MLE or MAP, depending on the 
assumed prior probabilities T’(±). 

We first simulate the nonadaptive decision rule. We 
choose times ranging from tf = 0.1ms to tf = 25 ms. 
For each trajectory, we use Eq. Q to calculate Aj^.. An 
error occurs if At / < Ath (Af^. > Ath) when the initial 
state is |-l-) (|—)). Here, Ath = P{~)/P{+) = is the 
optimal decision threshold. Averaging the errors over all 
trajectories for each state |±) gives the conditional error 
rate for each time tf. The error rate is then given 
by e = P(-|-)e+ -f P(—)e_. 

Next, we simulate the adaptive decision rule. We re¬ 
formulate the stopping condition in terms of the stopping 
probabilities p± = e'''±“'''“>/(l -I- Here, 


At/Ath 

1+At/Ath l + e^‘-N. 

is the probability of the initial state being |-|-) when the 
likelihood ratio is At = e^*. To vary the average read¬ 
out time, we choose stopping probabilities ranging from 
p+ = 0.9 to = 1 and from p_ = 0 to = 0.1. For 
each trajectory, we then use Eq. Q to find the first time 
t for which pt > p+, pt < P- or t > tM■ An error oc¬ 
curs if pt < 1/2 {pt > 1/2) when the initial state is |-|-) 
(|—)). Averaging the stopping times and errors over all 
trajectories for each state |±) gives the conditional aver¬ 
age readout time and the conditional error rate 
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FIG. 9. Evolution of the NV“ probability pt for a trajectory 
generated assuming |+) (solid blue) and |—) (dashed red) for 
the same parameters as in Fig. In this particular case, 
we assumed that P(+) = P{~)- The discontinuities in pt 
correspond to the detection of a photon. Data acquisition is 
stopped when pt first crosses the stopping thresholds p+ and 
P- (dotted black) or when t reaches the maximum acquisition 
time Im = 25 ms. 


for each pair of stopping probabilities {p+,p_}. The av¬ 
erage readout time and the error rate are calculated as 
T = P(-h)T(+)-hP(-)T(-) and e = P{+)e^+^+P{-)e^-\ 
respectively. We then numerically choose the pairs of 
stopping probabilities that minimize the error rate under 
the constraint of a constant average readout time. An 
example of the application of the adaptive decision rule 
is illustrated in Fig. for both initial states. 


Appendix D: Photon-counting method 


Here we summarize the analytical photon-counting 
readout model used in Ref. m- In this method, the 
number n of photons detected during the interval [0, tf] is 
compared to a threshold v. If n > v, we choose |-|-), while 
if n < we choose |—). To calculate the error rate, we 
need the conditional photon distributions P{n \±) . These 
can be obtained within our formalism via Eq. (B5). How¬ 
ever, from a computational point of view, it is more ef¬ 
ficient to use the expressions given in Ref. [T^]. These 
are 


tf. Each term in Eq. ( |D1[ ) is given by 
P{n,E\±) = V{n,fi[f)e-^^*f + 
r dtV{n, 

Jo tf t ( 02 ) 

P{n,0\±) = 

[ dtP{n, Jo( 2 a;t)e"’"±‘e"’"T(H-‘)^ 

Jo 

where Ik{z) is the fcth order modified Bessel function of 
the first kind 1391 and where we define 


= fJ—e 

nl 


xt = yr+r_t(t/ - 1). 
The conditional error rates are given by 
P{n\+), e(-) = Y 


n<.i/ 


n>h' 


Substituting Eq. (D2) into Eq. (D4), we obtain 

(±) — I 


where 




r dt q> {e, 

Jo tf -t 


d±) _ 
— 


(D3) 


(D4) 


(D5) 


(D6) 


Jo 

Here, the q> ’s are cumulative distribution functions of 
the Poisson distribution V{n,fj,): 


= Q{vP 1,^), 

= I - (5(i/-b 1,/r), 


(D7) 


where Q(a, z) is the regularized incomplete Gamma func¬ 
tion [53] . For a g iven time tf, we compute the in¬ 
tegrals in Eq. (|D6[ ) numerically to calculate the error 
rate e = P(-|-)e*TTq_ p(_)g(-) when P(-|-) is known, or 
e = -I- /2 when P(-l-) is unknown. We then 

increase the threshold starting from v = 0 until we find 
a local minimum in the error rate. We repeat the proce¬ 
dure at various times to plot the error rate as a function 
of time, as shown in Fig. 


P(n|±) = P(n,E|±)-bP(n,0|±). (Dl) 

Here, E (O) corresponds to an even (odd) number of 
ionizations or recombinations having occurred after time 


Appendix E: System rates and model verification 

In this section, we describe how the rates 7 ± and r± 
were extracted from experimental data and how our read¬ 
out model was validated. We first acquire 125 photon- 
count trajectories, each of duration 30 s and with time 
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bin <5t = 0.1 ms. We separate the trajectories into a cali¬ 
bration set (60 trajectories) and a testing set (65 trajec¬ 
tories). The rates are extracted using only the calibration 
set. The testing set is then used to independently verify 
the validity of the charge dynamics model extracted from 
the calibration set. 


1. Extraction of the system rates 


To determine the rate 7 ±, we split all trajectories of 
the calibration set into a total of 1800 subtrajectories of 
1 s, each of which is rebinned in bins of St = 10 ms. We 
then make a histogram of photon counts in all bins of 
10ms, as illustrated in Fig. [l0|(a). We fit the histogram 
to a mixture of Poisson distributions 


P{Sn) = V{Sn,"f+St)P{+) +'P{6n,^_St)P{-), (El) 

where T’(±) are prior probabilities for each charge state 
and where P{n,^) is given in Eq. (D3). The histogram 
is well fitted by Eq. (El), indicating that T±6t <C 1, i.e. 
r± ^ 100 Hz. From the means of the Poisson distribu¬ 
tions, we extract 7 + « 720 Hz and 7 _ « 50 Hz. When 
assuming no prior information on the state, the count 
threshold that best discriminates between NV“ and 
NV'^ in one time bin is given by 
or 


7-^ -7- 

'"(v) 


St. 


(E2) 


We find v « 2.5. Thus, we choose NV“ for 5n > 2 
and NV° for 6n < 2 , with a « 2 % rate of error for the 
simplified model of Eq. (El) (Fig. [^suggests that the 
actual error rate is closer to 3% at a measurement time 
of 10 ms). 

To determine F±, we use this thresholding procedure 
on the first 10 -ms bin of each 1 -s subtrajectory to post¬ 
select the initial state of the subtrajectory. This gives 
872 subtrajectories with the initial state identified as |-|-) 
and 2878 trajectories with the initial state identified as 
|—). We then obtain the average g±{t) of all trajectories 
for a given initial state |±), as illustrated in Fig. [To)(b). 
We fit both curves simultaneously [H] to th e rat e equa¬ 
tion model of the two-level fluctuator, Eq. (B2), which 
predicts 


0±{t) — Aj-e P 


(E3) 


where 


r = r. 


r_. 


^ r+(7_^t) + r_(7+Jt) 

T:+T_ ■ ^ ^ 


From the fitted values of F and B, we extract F_|_ « 3.6 Hz 
and F_ « 0.98 Hz. Note that in the future, the system 
rates could be extracted more accurately using machine 
learning |52j . The rates would then be chosen by training 
the MLE or MAP readout to minimize the empirical error 
rate. 




FIG. 10. (a) Histogram of photon counts per 10-ms time bin 
(black dots). The solid green curve is a fit to a mixture of 
Poisson distributions, Eq. (Ell. The dashed blue and dot¬ 
ted red curves are the Poisson distributions conditioned on 
states |-|-) and |—), respectively. We extract 7 + « 720 Hz and 
7 _ « 50 Hz from the means of the conditional distributions, 
(b) Average photon counts per 10 ms, Q±{t), for trajectories 
postselected on the initial states |-l-) (blue circles) and |—) (red 
squares). The solid black lines are a simultaneous ht of Q+{t) 
and g-{t) to Eq. (E3|. From the fit, we extract r+ « 3.6Hz 
and r_ « 0.98 Hz. The horizontal dotted black line is the 
steady-state average B obtained from the fit to Eq. (E3|. 


2. Verification of the readout model 

Finally, we verify the validity of the theoretical readout 
model in two complementary ways using the remaining 
testing set and the rates extracted from the calibration 
set. We first compare the experimental log-likelihood ra¬ 
tio distribution to the theoretical prediction. We then 
directly calculate the error rate by preparing and subse¬ 
quently measuring the charge state in postprocessing. 

To calculate the experimental probability density of 
the log-likelihood ratio, P{X), we split the testing set into 
78000 subtrajectories of duration 25 ms. Using Eq. 
we determine the value of the log-likelihood ratio A at the 
end of each subtrajectory. We use these values to calcu¬ 
late a histogram of P{X) using a bin size of SX = 0.1. 
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The resulting distribution is shown in Fig. ©a). We 
must compare the experimental distribution to the the¬ 
oretical prediction. We first calculate the expected con¬ 
ditional distribution Ptheo(A|±) using the Monte Carlo 
simulations described in Appendix]^ We then calculate 
the full distribution 


Ctheo(A) — ^ ^ Ctheo(A|?)Ftheo(^)- 


(E5) 


i=+.- 


We assume that the prior probabilities are given by their 
steady-state values, 


Ctheo( + ) — 1 ~ Ptheo{—) — 


r_ 


r+ + r_ 


(E 6 ) 


The theoretical distribution is compared to the experi¬ 
mental distribution in Fig. |ll[ a). The prediction repro¬ 
duces the experimental data without any fitting parame¬ 
ters. Our readout model can therefore be used to prepare 
the charge state with high fidelity in postselection. 

To directly obtain the experimental time dependence 
of the error rate, we split the testing set into 39000 sub¬ 
trajectories of duration 50 ms. We prepare the state by 
calculating the postmeasurement state p after the first 
25 ms of each subtrajectory, assuming a completely mixed 
initial state Pq = ( 1 / 2 , l/ 2 )'^(see the procedure described 
at the beginning of Appendix]^. If p+ > 1/2 (p+ < 1/2), 
we decide that the state is most likely |-|-) (|—)). In this 
way, we identify 8307 (30693) subtrajectories as being in 
the state |-|-) (|—)) after the first 25ms. We then use 
the MLE, adaptive MLE, and photon-counting methods 
to read out the state using the last 25 ms of each sub¬ 
trajectory. The outcome of the readout is compared to 
the preparation to give the measured conditional error 
rates e±. The measured error rate is then calculated as 
e = (e+ -I- e_)/2 and is plotted as a function of readout 
time in Fig.Jl^b). 

Figure [Tljbj shows that although the adaptive-decision 
speedup is unchanged, the measured error rate e is much 
larger than the theoretical prediction e. This could be the 
result of either imperfect preparation or imperfect mod¬ 
eling. We now show that this discrepancy is completely 
consistent with imperfect preparation. To account for 
preparation errors, we introduce an average preparation 
error rate rj, assuming equal prior probabilities. The true 
error rate is related to the measured error rate by the re¬ 
lation e = r]{l — e) -I- (1 — ri)e, or 


e = 


l-2ri' 


(E7) 


We simultaneously fit the experimental error-rate curves, 
using r] as the only fit parameter. This gives a best-fit 
value of 77 = 2.22%. We find a good agreement between 
theory and experiment, as shown in Fig. We conclude 
that the additional errors are likely due to preparation 
errors and not to imperfect modeling. We note that the 
transformation of Eq. (E71 does not affect the adaptive- 
decision speedup. 



Average readout time, T (ms) 

FIG. 11. (a) Experimentally measured probability density 

P(A) of the log-likelihood ratio A (black circles) for sub¬ 
trajectories of duration 25 ms. The theoretical prediction 
Ptheo(A) (dotted black) and the corresponding weighted con¬ 
ditional distributions Ptheo(A|-|-)Ptheo(+) (solid blue) and 
Ptheo(A|—)Ptheo(— ) (dashed red) are shown for comparison. 
The theoretical prior probabilities Ptheo(±) are taken to be 
the stationary values corresponding to the extracted rates, 
Eq. ( |E6[ ). (b) Experimental error rate e of the photon¬ 

counting method (blue circles), the nonadaptive MLE (orange 
squares), and the adaptive MLE (green triangles). The the¬ 
oretical error rate e of Fig. is reproduced for comparison 
(solid blue for the photon-counting method, dashed orange 
for the nonadaptive MLE, and dotted green for the adaptive 
MLE). The discrepancy between experiment and theory can 
be attributed to preparation errors [see Eq. (E71]. 
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